Correlated band theory of spin and orbital contributions to Dzyaloshinskii-Moriya 

interactions 
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A new approach for calculations of Dzyaloshinskii-Moriya interactions in molecules and crystals is 
proposed. It is based on the exact perturbation expansion of the total energy of weak ferromagnets 
in the canting angle with the only assumption of local Hubbard-type interactions. This scheme leads 
to a simple and transparent analytical expression for the Dzyaloshinskii-Moriya vector with a natural 
separation into spin and orbital contributions. The main problem was transferred to calculations 
of effective tight-binding parameters in the properly chosen basis including the spin-orbit coupling. 
Test calculations for La2Cu04 give the value of the canting angle in agreement with experimental 
data. 
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The Dzyaloshinskii-Moriya interactions (DMI) [1, 2] 
were introduced in the theory of weak ferromagnetism 
(WF) to explain the canting of the magnetic moments 
of some antiferromagnets (such as a-Fe203, M11CO3, 
C0CO3 and others [3(). It was shown later that the 
DMI are of crucial importance for many others classes 
of magnetic systems, such as spin glasses molecular 
magnets multiferroics syste ms (§1 , magnetic sur- 

faces and clusters on the surfaces [9l-tl3| , and Jahn- Teller 
In a sense, DMI is the simplest example 
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systems 

of relativistic magnetic interactions, since it appears al- 
ready in the first order in the spin-orbit coupling, whereas 
the magnetic anisotropy is at least of the second order Q ■ 
On the other hand, the DMI vanish for systems with in- 
version symmetry, which explains their special relevance 
for low symmetric cases such as molecules, clusters, sur- 
faces, and disorder systems. 

The microscopic origin of the DMI was clarified by 
Moriya for model systems For that he used the 

idea of Anderson[lJ] about the superexchange interac- 
tion mechanism. However, the original formulation of 
the Dzyaloshinskii-Moriya interaction is not suitable for 
quantitative calculations of the DMI-parameters for spe- 
cific compounds based on real electronic structures. Nu- 
merous attempts of more convenient a nd g eneral formu- 
lations have been made afterwards Til- 2lf. 

Yildirim et. al. 15J have developed a perturbative 



approach in the spin-orbit coupling for Mott insulators 
within the Hubbard-like model. A similar approach has 
been developed in Ref. [2l] within the LDA+U method. 
In Refs. mtH! the magnetic force theorem [22| has been 
applied, but only for spin rotations. The authors of 
Refl23l presented a computationally efficient method to 
determine the strength of the DMI from the spin-orbit 
induced corrections to the energy of long-ranged spin spi- 
rals. A general first-principle approach for the DMI was 



suggested in Ref. [l9| in a similar but fully relativistic 
formalism, that takes into account both spin and orbital 
contributions. This is probably the best possible way 
if one starts to calculate the DMI from the true non- 
collinear ground-state magnetic structure. 

The results of previous theoretical investigations (ljj 
have demonstrated that in the real transition metal com- 
pounds there are a lot of different microscopic mecha- 
nisms for the anisotropic exchange interactions. For in- 
stance, to take into account the metal-oxygen hybridiza- 
tion one should consider high-order hopping processes be- 
tween metal and oxygen orbitals. It strongly complicates 
the formulation and solution of the problem. In this pa- 
per we return to the original Anderson's idea [14] about 
the superexchange interaction in the Wannier function 
basis. We use the main advantage of such an approach 
which is that all the important hybridization effects can 
be captured by constructing the Wannier function. As we 
will show it simplifies dramatically the formalism without 
any essential loss of accuracy. 

Since the canting angles are normally quite small it 
allows us to proceed with the corresponding collinear 
structures and use advantage of first-order perturbation 
treatment for the magnetic torque. The application of 
the magnetic force theorem to equilibrium configurations, 
involves additional assumptions such as neglecting of ver- 
tex corrections [HI • First-order variation of the total en- 
ergy near the collinear states leads to an expression for 
the DMI formally exact in the many-body sense. 

We start with the general Hamiltonian of interacting 
electrons in a crystal: 

H = H t + H u 

= C 1^12 C 2 + \ C l C 2 U l23iC 3 C4, (1) 
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were 1 = (ii,mi,<ri) is the set of site (ix), orbital (mi) 
and spin (ui) quantum numbers and t\2 are hopping inte- 
grals that contain the spin-orbit coupling. These transfer 
couplings can be found by the Wannier-parameterization 
of the first- principle band structure with the spin-orbit 
coupling [25|. In this case the real space site-centered 
spinor Wannicr function can be written as 



W n (r) 



(2) 



where T is a lattice translation vector, t/^(r — T) are the 
site-centered spinor atomic-like orbitals (in our case they 
were linear muffin-tin orbitals (LMTO)) and w„ m t are 
expansion coefficients of the Wannier functions in terms 
of the corresponding LMTO orbitals. 

We will take into account only the local Hubbard-like 
interactions, keeping in H u only terms with i\ — ii = 
i% = «4. This assumption corresponds to the LDA+U 
Hamiltonian [261 ] that is also a starting point for the 
LDA+DMFT (Dynamical Mean-Field Theory) [27l.l28jj. 
It is crucially important for the later consideration that 
the interaction term H u is supposed to be rotationally 
invariant. 

Let us start with a collinear magnetic configuration 
(e.g. the Neel antiferromagnetic state), which is close to 
the real ground state (weak ferromagnet), but does not 
coincide with it due to the Dzyaloshinskii-Moriya inter- 
actions. The phenomenological Hamiltonian of the DMI 
is given by 



DM 



e 7 ; X €j 



(3) 



where is a unit vector in the direction of the i-th site 
magnetic moment and is the Dzyaloshinskii-Moriya 
vector. We analyze the magnetic configuration that is 
slightly deviated from the collinear state, 



Si = i]ie + [Sifi x rfceb], 



(4) 



where rji = ±1, e*o is the unit vector along the vector 
of antiferromagnetism, and 50i are the vectors of small 
angular rotations. 

Substituting Eq. (|4]) into Eq. ([3]) one finds for the 
variation of the magnetic energy: 



5E = y^AjQftpi - 80 j)- 



(5) 



on the total energy; here J; = Li + Si is the total moment 

operator, Li and Si are the orbital and spin moments, re- 
spectively. We would like to stress that the operator Ri 
acts on ith Wannier state. In the Supplementary mate- 
rials [29[ we demonstrate that the rotation of the orbital 

part of Ji in Wannier function basis results in indepen- 
dent rotations of the atomic orbital moments. 

The interaction part of the Hamiltonian H u is rota- 
tionally invariant and is not changed under this transfor- 
mation, opposite to the hopping part H t ■ 

SH t = J2 4( S ^ftij + iijSRj)^ 
= -i^^cf(80iJiiij - iijJjS0j)cj 

ij 



(7) 



Assuming that Ji 
takes the form 



Jj = J the change of the total energy 



5E = - - ^2(5<pi - 5yj)Tr m , a {cj[J,tij] 



"o ^2( S <Pi + S <Pj) Tr m,A c t[Jiiij]-Cj), (8) 



where Tr m a is a trace over orbital (m) and spin (a) 
quantum numbers. 

The first term in the right-hand side of Eq. ([SJ is 
responsible for relative deviations of the magnetic mo- 
ments on sites i and j (DMI) whereas the second one is 
related with the rotation of the magnetic axis as a whole 
(magnetic anisotropy). Assuming that S0i = 50 is inde- 
pendent on site index one finds the following expression 
for the magnetic anisotropy torque 



5E 

Sep 



-i ^ Tr m . a (c+ [J, tij]- Cj ) . 



(9) 



In contrast with the previous results [17J, UM the expres- 
sion ^) contains both spin and orbital contributions. Ap- 
plication of this expression to real systems will be con- 
sidered elsewhere. Here we will focus on the DMI. 
Comparing Eq. (j8]) with Eq. ([5]) one finds 



Now we should calculate the same variation for the 
microscopic Hamiltonian (JlJ. Similar to the procedure 
used in Ref. [2J to derive exchange interactions for the 
LDA+DMFT approach, we consider the effect of the lo- 
cal rotations 



JSipiJi 



(6) 



Di 



'[Jiiijj+Cj) 



-Trm^NjiiJJijU^O) 



where Nji — (cfcj) = — i J^ImGji(E)dE is the oc- 
cupation matrix and G is the Green function, Ef is 
the Fermi energy. Assuming that the occupation ma- 
trix is known exactly the expression for DMI (ITU1) is ex- 
act due to the Hellmann-Feynman theorem. Note that 
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the occupation matrix is calculated in the correspond- 
ing collinear states, which strictly speaking can be done 
self-consistently only within constrained calculations [3(| ■ 

Using the decomposition of the total moment J into or- 
bital and spin moments, we have a natural representation 
of the Dzyaloshinskii-Moriya vector (ITUl) as a sum of the 
orbital and spin contributions which are related with the 
rotations in orbital and spin space, respectively. 

To test the developed method we consider weak fer- 
romagnetism phenomena which result from DMI. The 
problem of the theoretical description of weak ferromag- 
netism in antiferromagnets can be solved by calculating 
the canting angle. As an example of the system demon- 
strating weak ferromagnetism we have chosen La2CuC>4 
in the low-temperature orthorhombic phase presented in 
Fig.l. For this system the Wannier functions can be qual- 
itatively analyzed by using a one-band Hubbard model 
with the spin-orbit coupling proposed in Refs. 



netic moments were fixed along z-axis. The computa- 
tional details are the same as in Ref. [U The obtained 
results are presented in Table I. 
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H = ^ C ta( tS af3 +i\ijV a p)Cjp +U }^ Tt^n^. (11) 



Here t is a nearest-neighbor hopping parameter and the 
vector Xij depends on the tilting pattern of oxygen oc- 
tahedra surrounding the copper atom. Substituting Eq. 
([TTjl into Eq. (JTOj) we obtain 



FIG. 1: Magnetic structures of La2Cu04. Black and grey 
arrows denote non-collinear and fixed collinear ground 
states, respectively. Grey and yellow circles represent 
oxygen atoms which are coming out of and going into 
the copper-oxide plane. 

In order to calculate the Dzyaloshinskii-Moriya inter- 
action (|10[) one needs to define the occupation matrix Nji 
and the kinetic part iy of the Hamiltonian in a Wan- 
nier function basis. The construction of a reliable Wan- 



Dij = XijTr a Nji. 



(12) 



nier basis can be performed in different schemes [36H39| . 
As a first attempt we will use the simplest choice re- 
lated with orthogonalized minimal basis LMTO scheme 
36l 38 1 including the spin-orbital coupling. Note that, 
since the trace in Eq. (10) is only over orbital indices 



Therefore the symmetry of the Dzyaloshinskii-Moriya and no * on the site „° nes : J he resulting Aj is Wannier- 
vector is fully described by the vector Ay. If there is 
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the inversion center between copper atoms than Ay 
and DMI vanishes. 

Using the definitions introduced in Refs 
obtain A12 = A i4 = i(Ai -A 2 , Ai + A 2 , 0) and A13 = A15 = 
|(Ai — A2, — Ai — A2, 0). It is easy to show that the total 
Dzyaloshinskii-Moriya vector Y^j=2 D\j has only nonzero 
x component. It means that the canting exists if the 
magnetic moments are in yz plane. This fully agrees 
with the results of previous works [31, 33, |34|. 



Unfortunately, the performed microscopic analysis is 
only qualitative. The most straightforward way to obtain 
a reliable numerical estimation of the canting angle would 
be to perform relativistic first-principle calculations for 
the corresponding noncollinear magnetic structure. How- 
ever, this is a very challenging computational problem. 
Since in the case of 3d-metal compounds the spin-orbit 
coupling is small one should run many thousands of it- 
erations to relax the magnetic structure [35| . On the 
other hand, to solve the weak ferromagnetism problem 
we only need to know the magnitude and the direction of 
the Dzyaloshinskii-Moriya vector. It can be done by us- 
ing the developed method for the fixed collinear magnetic 
configuration. Such an approach seems to be preferable 
since it requires much less computational efforts. 

We have performed the LDA+U+SO calculations for 
the collinear antiferromagnetic structure where the mag- 



gauge-dependent. The Wannier functions is normally 
constructed using a truncated basis, therefore it will be 
important to investigate in the future a sensitivity of the 
results with respect to a choice of the Wannier states. 

Since in our investigation we used the atomic sphere 
approximation (LMTO-ASA)[36] it was natural to asso- 
ciate Nji and £y with the occupation matrix and kinetic 
energy of the 3d states of the copper atoms. We consider 
such an approximation as a reasonable one since the mag- 
netic moments in La2Cu04 are due to the 3d states of 
copper [i^]. Another fact supported our approximation 
is good agreement between the isotropic exchange inter- 
action calculated by using the Green's function method 
in the framework of LMTO- ASA and the model kinetic 

--ii Ell 



exchange estimated as J,. 



TABLE I: Calculated magnitude (in (ib) and orientation of 
the spin and orbital copper moments in La2Cu0 4 . 



Atom Spin moment Orbital moment 

1 0.65 x (0, 0, -1) 0.04 x (0, 0, -1) 

2 0.65 x (0, 0, 1) 0.04 x (0, 0, 1) 



The DMI parameters between neighboring copper 
atoms calculated via Eq. (fT0| are presented in Ta- 
ble II. One can see that the orbital contribution to the 
Dzyaloshinskii-Moriya interaction is one order of magni- 
tude larger than the spin one. The obtained magnetic 
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TABLE II: Different contributions to the Dzyaloshinskii- 
Moriya vectors (in meV). 



Ri 



(1.2) (-0.005;-0.006; 0) (-0.07; -0.03; 0) 

(1.3) (-0.005; 0.006; 0) (-0.07;0.03; 0) 

(1.4) (-0.005; -0.006; 0) (-0.07;-0.03; 0) 

(1.5) (-0.005; 0.006; 0) (-0.07;0.03; 0) 



torque is directed along x axis. This agrees with the re- 
sults of our microscopic analysis. Summarizing all the 
vectors we can calculate the canting angle of the mag- 
netic moment which is given by 
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0.005, 



(13) 



where the total exchange interaction Y^j Jij was taken 
to be 58.3 meV [21] . The obtained value of the canting 
angle is in a reasonable agreement with the experimental 
estimate of 0.003 0. 

To conclude, we have proposed a new method for 
calculation of the Dzyaloshinskii-Moriya interaction pa- 
rameters which, conceptually, is much simpler than ap- 
proaches known before. This method represents, in a 
natural way, the Dzyaloshinskii-Moriya vector as a sum 
of the spin and orbital contributions which may give a 
deeper insight into microscopic mechanisms of the DMI 
for a given system. In our approach, the crucial point is 
the construction of a reliable tight-binding parameteri- 
zation of the Hamiltonian with the spin-orbit interaction 
taken into account. We have performed the correspond- 
ing calculations for the weak ferromagnet La2Cu04, and 
the results look quite promising. 
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